Living in a Void: Testing the Copernican Principle with Distant Supernovae 
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A fundamental presupposition of modern cosmology is the Copernican Principle; that we are not 
in a central, or otherwise special region of the Universe. Studies of Type la supernovae, together 
with the Copernican Principle, have led to the inference that the Universe is accelerating in its 
expansion. The usual explanation for this is that there must exist a 'Dark Energy', to drive the 
acceleration. Alternatively, it could be the case that the Copernican Principle is invalid, and that 
the data has been interpreted within an inappropriate theoretical frame-work. If we were to live in 
a special place in the Universe, near the centre of a void where the local matter density is low, then 
the supernovae observations could be accounted for without the addition of dark energy. We show 
that the local redshift dependence of the luminosity distance can be used as a clear discriminant 
between these two paradigms. Future surveys of Type la supernovae that focus on a redshift range 
of ~ 0.1 — 0.4 will be ideally suited to test this hypothesis, and hence to observationally determine 
the validity of the Copernican Principle on new scales, as well as probing the degree to which dark 
energy must be considered a necessary ingredient in the Universe. 



The concordance model of the Universe combines two 
fundamental assumptions. The first is that space-time 
is dynamical, obeying Einstein's Equations. The second 
is the 'Cosmological Principle', that the Universe is then 
homogeneous and isotropic on large scales - a general- 
isation of the Copernican Principle that "i/ie Earth is 
not in a central, specially favored position" [jj. As a re- 
sult of these two assumptions we can use the Freidmann- 
Robertson- Walker (FRW) metric to describe the geome- 
try of the Universe in terms of a single function, the scale 
factor a{t), which obeys 



H' = 



8ttG 



k 



(1) 



where H = a/ a is the Hubble rate, p is the energy den- 
sity, k is the (constant) curvature of space, and overdots 
denote time derivatives. The scale factor can then be 
determined by observing the 'luminosity distance' of as- 
trophysical objects. At small z = ao/a{t) — 1 this is given 
by 



HqDl ~ cz + -(1 - qo)cz^, 



(2) 



where q = — aa/d^ is the deceleration rate, and subscript 
denotes the value of a quantity today. Recent mea- 
surements of (z, Dl) using high redshift. Type la Super- 
novae (SNe) have indicated that an < 0, i.e. the Universe 

i1 rTi 
is accelerating in its expansion [2, |3[. Accelerating ex- 
pansion is possible in an FRW universe if a fraction of 
p is in the form of a smoothly distributed and gravita- 
tionally repulsive exotic substance, often referred to as 
Dark Energy [J| . The existence of such an unusual sub- 
stance is unexpected, and requires previously unimagined 
amounts of fine-tuning in order to reproduce the obser- 
vations. Nonetheless, dark energy has been incorporated 
into the standard cosmological model, known as ACDM. 
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An alternative to admitting the existence of dark en- 
ergy is to review the postulates that necessitate its intro- 
duction. In particular, it has been proposed that the SNe 
observations could be accounted for without dark energy 
if our local environment were emptier than the surround- 
ing Universe, i.e. if we were to live in a void (a, la, |7|. 
This explanation for the apparent acceleration does not 
invoke any exotic substances, extra dimensions, or modi- 
fications to gravity ~ but it does require a rejection of the 
Copernican Principle. We would be required to live near 
the centre of a spherically symmetric under-density, on 
a scale of the same order of magnitude as the observable 
Universe. Such a situation would have profound conse- 
quences for the interpretation of all cosmological obser- 
vations, and would ultimately mean that we could not 
infer the properties of the Universe at large from what 
we observe locally. 

Within the standard inflationary cosmological model 
the probability of large, deep voids occurring is extremely 
small. However, it can be argued that the centre of a 
large underdensity is the most likely place for observers 
to find themselves Q. In this case, finding ourselves in 
the centre of a giant void would violate the Copernican 
principle, that we arc not in a special place, but it may 
not violate the Principle of Mediocrity, that we are a 'typ- 
ical' set of observers. Regardless of what we consider the 
a priori likelihood of such structures to be, we find that 
it should be possible for observers at their centre to be 
able to observationally distinguish themselves from their 
counterparts in FRW universes. Living in a void leads to 
a distinctive observational signature that, while broadly 
similar to ACDM, differs qualitatively in its details. This 
gives us a simple test of a fundamental principle of mod- 
ern cosmology, as well as allowing us to subject a possible 
explanation for the observed acceleration to experimental 
scrutiny. 

Some efforts have gone into identifying the observa- 
tional signatures that could result from living in a void. 
The cosmic microwave background (CMB) supplies us 
with the tight constraint that we must be within 15 Mpc 
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FIG. 1: Some representative forms for the curvature function, fc(r), and the corresponding distance moduli, normalised to an 
empty Milne cosmology, Adm. Void (a) is a Gaussian in r (with FWHM ro), void (b) is a fc oc exp{— c|r|'^}, and void (c) is 
fc oc {1 — I tanh(r)|}, all normalised to the curvature minimum ko. The three diagrams on top indicate the spatial distribution 
of energy density in the three voids at some sample time (solid line for void (a), dashed line for void (b) and dotted line for 
void (c)). The ascending thick solid line on the right hand plot corresponds to a universe containing dark energy only (de Sitter 
space), and the descending one to a flat universe containing only dust (an Einstein-de Sitter universe). 



of the center of the void [9| . There have also been some 
attempts at calculating predictions for CMB anisotropics 
and large scale-structure [lO, [lllll3 , as well as the kine- 
matic Sunyaev-Zeldovich effect [ij] . 

General Relativity allows a simple description of time- 
dependent, spherical symmetric universes: the Lemaitre- 
Tolman-Bondi (LTB) models [11 [H, [H, whose hne- 
element is 



ds^ = -dt^ 



a\{t,r)dr'^ 
1 — k{r)r^ 



+ al(t,rydQ^ 



(3) 



where 02 = (rai)', and primes denote r derivatives . The 
old FRW scale factor, a, has now been replaced by two 
new scale factors, oi and 02, describing expansion in the 
directions tangential and normal to the surfaces of spher- 
ical symmetry. These new scale factors are functions of 
time, t, and distance, r, from the centre of symmetry, 
and obey a generalization of the usual Friedman equa- 
tion, ([1]), such that 



SttG^ 



fc(r) 



(4) 



Here p = ra{r)/a\, and is related to the physical energy 
density by p = p-\-p'rai/'ia2- The two free functions, k{r) 
and ■m{r), correspond to the curvature of space, and the 
distribution of gravitating mass in that space. We choose 
initial conditions such that the curvature is asymptoti- 
cally flat with a negative perturbation near the origin, 
and so that the gravitational mass is evenly distributed. 



As the space-time evolves the energy density in the vicin- 
ity of the curvature perturbation is then dispersed, and 
a void forms. Observations of distant astro-physical ob- 
jects in this space-time obey a distance-redshift relation 



Dj 



(1 + z)'^rEai{tE,rE) 



where 



1 + z = cxp 



{airy 



zdr 



(5) 



(6) 



and subscript E denotes the value of a quantity at the 
moment the observed photon was emitted. This expres- 
sion is modified from equation ^, allowing for the pos- 
sibility of apparent acceleration without dark energy. 

We find that the form of the void's curvature profile 
is of great importance for the observations made by as- 
tronomers at its centre. In Figure [1] we plot some simple 
curvature profiles, together with the corresponding dis- 
tance moduli as functions of redshift (distance modulus, 
Adm, is defined as the observable magnitude of an as- 
trophysical object, minus the magnitude such an object 
would have at the same redshift in an empty, homoge- 
neous Milne universe). It is clear from Figure [1] that 
for the void models there is a strong correlation between 
k{r) and Adm; at low redshifts Adm(z) traces the shape 
of k{r). Hence, for a generic, smooth void Adm starts 
off with near zero slope, where it is locally very similar 
to a Milne universe, it then increases at intermediate z, 
and later drops off like an Einstein-de Sitter universe. 



For ACDM, we have Adm 



'jQqz at low z, i.e. 




non-zero slope. Thus, although one can always find a 
void profile that will mimic ACDM ^7|] , such a void will 
have a curvature profile that is strongly cusped, and non- 
differentiable, at z = [l8| (see the dotted line in Figure 
fflor Lis])- Conversely, any generalized dark energy model 
capable of producing a flat Adm at low z would be re- 
quired to change the equation of state extremely rapidly 
between z ~ 0.5 and 0.1. We therefore have a definitive 
way to distinguish between a realistic smooth void model, 
and ACDM. 
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FIG. 2: The current best fit ACDM and Gaussian void mod- 
els as dashed and solid lines, with triangular and square 
data points, respectively. The data shown here is a compi- 
lation of 115 low and high-2 SNe from the SNLS, fitted using 
SALT [20|]. For illustrative purposes we have binned the re- 
sults with 10 SNe per bin (except the last one, which contains 
5). Due to the uncertainty in the 'nuisance' parameters of the 
calibration the data points and the error-bars shift when fit- 
ting for different models. 

We will now compare the smooth void model to the 
first-year SNe Legacy Survey (SNLS) data, consisting 
of 115 SNe ^ cahbrated with the SALT light-curve 
fitter [m, and contrast it with ACDM. We use the 
Bayesian information criterion as a figure of merit (see 
e.g. [22, 1231), and assume one model to be decisively 
favored over the other if 



\A\nE\K |A(ln/:„ 



P 



lniV)| > 5, 



(7) 



where E is the evidence for a model, Cmax is the max- 
imum likelihood of a model, given a data set, p is the 
number of parameters in the model, and N is the number 
of data points. This criterion corresponds to one model 
being ~150 times more likely than the other. The mini- 
mal void model under consideration has 6 parameters: 2 
to parametrize a Gaussian fc, and 4 'nuisance' parameters 
required to calibrate the SNe data. These are absolute 
magnitude, intrinsic error, and the colour and stretch pa- 
rameters used in light curve fitting, {Mg, crint, a, /3}. As- 
suming spatial flatness, ACDM requires 5 parameters: 1 
specifying the fraction of dark energy, and the same 4 nui- 
sance parameters. In a more comprehensive study it may 
be preferable to perform a full Bayesian evidence analy- 
sis, with suitable priors [23|. In the interests of brevity, 
and to avoid a lengthy discussion of prior probabilities, 
we have refrained from this for now. 

In Figure [2] we show the SNLS data with the two best 
fit models. Both have similar goodness of fit, but one can 
discern a qualitative difference between them, which will 



be distinguishable with future surveys. The best fitting 
void is 71 ± 7% underdense at its centre, and has a scale 
corresponding to 850 ± 170/i~^Mpc today. This is of the 
order expected to produce a feature in Adm on a scale 
of z ~ 0.6, and large enough to avoid strong constraints 
from galaxy surveys that extend to z ~ 0.1. On the 
other hand, the best fitting ACDM model contains 74 ± 
4% dark energy, and fits the data slightly better with 
|Alni?| ~ 2.7. Thus, while the current data marginally 
prefers ACDM, it is not yet able to distinguish between 
the two models decisively. 

One will, of course, be interested in results from other 
SNe compilations. Using the Riess gold data [2J, with 
the MLCS2k2 light-curve fitter, we find our basic results 
do not change significantly, with a ACDM model still 
being marginally preferred. Thus our analysis does not 
appear to be substantially effected by the apparent sys- 
tematic error that led to the identification of the 'Hubble 
Bubble' anomaly [2^. The 'Union' data of [2y| is the 
largest compilation of SNe fitted for with the more con- 
servative SALT fitter, and for these 315 SNe (including 
the 'outliers') we find the void model is marginally pre- 
ferred over ACDM, with |Alni5| ~ 2.5 in favor of a void. 
It therefore appears neccessary to obtain more data, in 
order to be able to decisively distinguish the two models. 
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FIG. 3: The best fit ACDM and Gaussian void models as 
dashed and solid lines, with triangular and square data points, 
respectively, for an example of 700 SNe simulated from a 
ACDM model using the SALT light curve fitter. The shape of 
the redshift distribution used here is similar to that expected 
from the 2000 JDEM SNe, with an extra 300 at low z. We find 
that 700 SNe, with this redshift distribution, can decisively 
recover the ACDM model 99% of the time, as it becomes 
evident that a void model cannot mimic the low redshift be- 
haviour of a ACDM cosmology. For illustrative purposes, we 
have binned in groups of 60. 

Due to the different strategies and technologies, SNe 
surveys typically target either low or high redshifts. This 
has lead to a dearth of SNe at0.1<z<0.4- exactly 
the location where there is the greatest qualitative differ- 
ence between the two models. Future SNe surveys, with 
a redshift coverage in this region, will do better. As an 
example of the future constraints we can expect to gain, 
we consider the JDEM missions which expect to observe 
-2000 high redshift SNe in the interval 0.1 < z < 1.7, 
with an expected smooth distribution [27J. At very low 
redshifts (0.03 < z < 0.1) a further -300 SNe can be 
expected to be observed by other projects. We consider 
data simulated from a ACDM model, using the SALT fit- 



ter, and a Gaussian void model with the same parameter 
values as before. In both cases we find that these 2300 
SNe are sufficient to decisively recover the correct under- 
lying model. For the void model we find A\nE = 89± 12, 
while for the ACDM model AlnE ^ 46± 10. These are 
considerably better than the decisive benchmark, which 
was satisfied by all 1000 of our simulations (except one in 
the ACDM case). We can therefore say with confidence 
that the upcoming SNe data will be able to distinguish 
ACDM from a void model, independent of which is re- 
sponsible for the apparent acceleration. 

Given that --2300 SNc from a JDEM mission will do a 
superfluous job, we will now estimate the minimum num- 
ber of SNe required to meet the bound in Equation ([7]). 
Using the same redshift distribution as before [27|) we 
now consider different numbers of SNe. In the case of 
the void cosmology we find that with ^170 SNe 50% of 
our 1000 simulations recovered the void model decisively, 
while with ~480 SNe 99% of the simulations could do so. 
Similarly, in the case of a ACDM cosmology we found 
that with '--^ISO and ~700 SNe we could decisively recover 
the correct model in 50% and 99% of the simulations, re- 
spectively. This is illustrated in Fig. 3. These projections 
are much lower than the number of SNe expected to be 
observed by the next generation of SNe surveys, and we 
should therefore expect less ambitious projects to able to 
distinguish between the two models. In fact, this may be 
possible soon, with data from the Sloan Survey and third 
year SNLS expected to be released imminently. The re- 
sults obtained here depend on the details of the future 
surveys. Tailoring these to the specific regions where the 
two models differ most would undoubtedly give decisive 
results sooner. 

Consider now the effect of varying the redshift distri- 
bution, rather than the overall number. Adding a further 



300 SNe to the SNLS data, and letting them have a Gaus- 
sian redshift distribution with a^ — 0.1 and a mean that 
can vary, we simulate the data from void and ACDM 
models. We find that in a void Universe the optimal 
place to search for these extra 300 SNe is at z ^ 0.3 (lim- 
iting ourselves to a mean redshift of less than 1), where 
the average Aln£max '^-^ 5. In a ACDM universe the SNe 
would be better placed a little lower, at z '^ 0.1, and in 
this case Aln£niax ~ 2.5. 

We emphasise that similar results to those presented 
above should be obtainable for any smooth void model. 
For example, repeating our analysis for void (b), in Fig. 
[TJ gives almost identical results to void (a). The rea- 
son for this is that all smooth voids display the same 
qualitative behaviour of having a fiat Adm(z) profile at 
low-z. It is this qualitative difference that allows the void 
models to be so easily distinguished from ACDM, and as 
all smooth voids display this low-z behaviour, we expect 
the results we have presented to be broadly generalizable 
to all voids. Of course, it may be possible to imagine 
anomolous cases. This will be considered further in a 
more extended future publication. 

Two very different paradigms have been invoked to ex- 
plain the current observation of an apparently accelerat- 
ing Universe, depending on whether we invoke or reject 
the Copernican Principle. We have shown that in the 
coming years it will be possible to experimentally distin- 
guish between these two scenarios, allowing us to exper- 
imentally test the Copernican Principle |2^, [2^, [3(J| , as 
well as determine the extent to which Dark Energy must 
be considered a neccessary ingredient in the Universe. 
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